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Abstract 

Qh' We present a formulation of the determination of the impedance 

r~| ■ between any two nodes in an impedance network. An impedance 

network is described by its Laplacian matrix L which has generally 
complex matrix elements. We show that by solving the equation 
Lu a = \ a u* a with orthonormal vectors u a , the effective impedance 
between nodes p and q of the network is Z pq = J2a( u ap ~ Uaq) 2 /^a 
^ ' where the summation is over all X a not identically equal to zero and 

■ u ap is the p-th component of u a . For networks consisting of induc- 

tances L and capacitances C, the formulation leads to the occurrence 
of resonances at frequencies associated with the vanishing of X a . This 
curious result suggests the possibility of practical applications to res- 
onant circuits. Our formulation is illustrated by explicit examples. 
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1 Introduction 



A classic problem in electric circuit theory that has attracted attention 
from Kirchhoff 's time to the present is the consideration of network 
resistances and impedances. While the evaluation of resistances and 
impedances can in principle be carried out for any given network using 
traditional, but often tedious, analysis such as the Kirchhoff's laws, 
there has been no conceptually simple solution. Indeed, the problem 
of computing the effective resistance between two arbitrary nodes in 
a resistor network has been studied by numerous authors (for a list 
of relevant references on resistor networks up to 2000 see, e.g., J2J). 
Particularly, an elementary exposition of the material can be found in 
Doyle and Snell 

However, past efforts prior to 2004 have been focused mainly on 
regular lattices and the use of the Green's function technique, for 
which the analysis is most conveniently carried out when the network 
size is infinite [21 II]- Little attention has been paid to finite networks, 
even though the latter are those occurring in applications. Further- 
more, there has been very few studies on impedance networks. To be 
sure, studies have been carried out on electric and optical properties of 
random impedance networks in binary composite media (for a review 
see [3]) and in dielectric resonances occurring in clusters embedded in 
a regular lattice 0. But these are mostly approximate treatments on 
random media. More recently Asad et al. [7j evaluated the two-point 
capacitance in an infinite network of identical capacitances. When all 
impedances in a network are identical, however, the Green's function 
technique used and the results are essentially the same as those of 
identical resistors. 

In 2004 one of us proposed |Hj a new formulation of resistor net- 
works which leads to an expression of the effective resistance between 
any two nodes in a network in terms of the eigenvalues and eigenvec- 
tors of the Laplacian matrix. Using this formulation one computes the 
effective resistance between two arbitrary nodes in any network which 
can be either finite or infinite This is a fundamentally new formu- 
lation. But the analysis presented in |B] makes use of the fact that, 
for resistors the Laplacian matrix has real matrix elements. Conse- 
quently the method does not extend to impedances whose Laplacian 
matrix elements are generally complex (see e.g., [0]). In this paper we 
resolve this difficulty and extend the formulation of |H] to impedance 
networks. 
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Consider an impedance network C consisting of N nodes numbered 
a = 1, 2, M. Let the impedance connecting nodes a and (3 be 



z a /3 = zp a = r a p + i x a(3 



(1) 



where r a p = rp a > is the resistive part and x a p = xp a the reactive 
part, which is positive for inductances and negative for capacitances. 
Here, i = often denoted by j = y—1 in alternating current (AC) 
circuit theory In this paper we shall use i and j interchange- 
ably. The admittance y connecting two nodes is the reciprocal of the 
impedance. For example, y a p = yp a = l/z a p. 

Denote the electric potential at node a by V a and the net current 
flowing into the network (from the outside world) at node a by I a . 
Both V a and I a are generally complex in the phasor notation used in 
AC circuit theory [Sj. Since there is neither source nor sink of currents, 
one has the conservation rule 



a=l 



0. 



The Kirchhoff equation for the network reads 



where 



/ 2/1 



LV = I 



"2/12 



(2) 



(3) 



with 



-2/21 V2 

V -ym -vm 



y a = 2/ q /9' 

(3=1 {p^a) 



-yw \ 

-y2jv 

yx J 



(4) 



(5) 



is the Laplacian matrix associated with the network C In Q, V and 
/ are N- vectors whose components are respectively V a and I a . 

Here, we need to solve (J3J) for V for a given current configuration 
/. The effective impedance between nodes p and q, the quantity we 
wish to compute, is by definition the ratio 



V„ - V„ 



(6) 



3 



where V p and V q are solved from © with 

la = I {S ap -S aq ). (7) 

The crux of matter is to solve the Kirchhoff equation (J3J) for I 
given by (J7J). The difficulty lies in the fact that, since the matrix L is 
singular, the equation Q cannot be formally inverted. 

To circumvent this difficulty we proceed as in jS] to consider instead 
the equation 

L(e)V(e) = I (8) 

where 

L(e)=L + eI (9) 

and I is the identity matrix. The matrix L(e) now has an inverse and 
we can proceed by applying the arsenal of linear algebra. We take the 
e — > limit at the end and do not expect any problem since we know 
there is a physical solution. 

The crucial step is the computation of the inverse matrix L _1 (e). 
For this purpose it is useful to first recall the approach for resistor 
networks. 

In the case of resistor networks the matrix L(e) is real symmet- 
ric and hence it has orthonormal eigenvectors ip a (e) with eigenvalues 
= A Q + e determined from the eigenvalue equation 

L(e)Va(e) = Aa(e)Va(e), i = 1,2, ...,J\f. (10) 

Now a real hermitian matrix L(e) is diagonalized by the unitary trans- 
formation TJT(e) L(e) U(e) = A(e), where U(e) is a unitary matrix 
whose columns are the orthonormal eigenvectors ip a (^) and A(e) is a 
diagonal matrix with diagonal elements A Q (e) = X a + e. The inverse 
of this relation leads to L _1 (e) = U(e) A _1 (e) U^(e). 1 In this way we 
find the effective resistance between nodes p and q to be jS] 

M 1 2 

R pq = T^ap - ^«g| 2 (11) 

where the summation is over all nonzero eigenvalues, and ip ap is the 
pth component of t/j a (0). Here the a = 1 term in the summation with 

1 The equivalent of the method we use in obtaining is known in mathematics 
literature as the pseudo-inverse method (see, e.g., 
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Ai(e) = e and ipi p (e) = 1/y/Af drops out (before taking the e — > 
limit) due to the conservation rule (|2"|). It can be shown that there 
is no other zero eigenvalue if the network is singly-connected. The 
relation is the main result of [S|. 

2 Impedance networks 

For impedance networks the Laplacian matrix L is symmetric and 
generally complex and thus 

Lt = L* + L 

where * denotes the complex conjugation and f denotes the hermitian 
conjugation. Therefore L is not hermitian and cannot be diagonalized 
as described in the preceding section. 

However, the matrix L is always hermitian and has nonnegative 
eigenvalues. Write the eigenvalue equation as 

L ip a = a a ip a , cr Q > 0, a = 1,2, ...,J\f. (12) 

One verifies that one eigenvalue is o\ = withal = {1, 1, 1} T /vW, 
where the superscript T denotes the transpose. For complex L there 
can exist other zero eigenvalues (see below). 

To facilitate considerations, we again introduce L(e) as in and 
rewrite ()12j) as 

Lt(e)L(e)^ a (e) = a a (e)t/j a (e) , a a {e) > 0, a = l,2,..,AA, (13) 

where e is small. Now one eigenvalue is <Ji(e) = e 2 with V'i( e ) = 
{1, 1, l} T /\/7\f. For other eigenvectors we make use of the theo- 
rem established in the next section (see also |12| ) that there exist M 
orthonormal vectors u a {e) satisfying the equation 

L(e)« a (e) = A a (e)«*(e), a = l,2,...,M (14) 

where 

Aa(e) = ^{e)e w ^\ 9 a (e) = real. (15) 
Particularly, we can take 

A 1 (e) = y^)=€ ) 9 1 (e) = 0. (16) 
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Equation (|14[) plays the role of the eigenvalue equation 1)10(1 for resis- 
tors. 

We next construct a unitary matrix U(e) whose columns are u a {e). 
Using (|14j) and the fact that L(e) is symmetric, one verifies that L(e) 
is diagonalized by the transformation 

U T (6)L(6)U(6)=A(6) 

where A(e) is a diagonal matrix with diagonal elements A a (e). The 
inverse of this relation leads to 

L~ 1 (e) = \J(e)A- 1 (e)V T (e). (17) 

where A~ 1 (e) is a diagonal matrix with diagonal elements l/A Q (e). 
We can now use (fT7|) to solve (jHJ) to obtain, after using ©, 




where u ap is the pth component of the orthonormal vector u a (e). 

Now the term a = 1 in the summation drops out before taking the 
limit just like in the case of resistors^ since Ai(e) = e and u\ p {e) = 
ui q (e) = constant. If there exist other eigenvalues A Q (e) = e with 
Uap(c) 7^ constant, a situation which can happen when there are pure 
reactances L and C, the corresponding terms in (|18|) diverge in the 
e — > limit at specific frequencies u in an AC circuit. Then one 
obtains the effective impedance 

N 1 2 

Z pq = T-( u ap ~ u aq ) , if X a / 0, a > 2 

a=2 Aa V 7 

= oo, if there exists A Q = 0, a > 2. (19) 

Here u ap = u Qp (0). The physical interpretation of Z = oo is the 
occurrence of a T6son(iTicc in an AC circuit at frequencies where l a — 
0, meaning it requires essentially a zero input current I to maintain 
potential differences at these frequencies. 

The expression ()19|) is our main result for impedance networks. 

In the case of pure resistors, the Laplacian L(e) and the eigenvalues 
l a (e) in (|1U|) are real, so without the loss of generality we can take ip a (e) 
to be real (see Example 3 in Sec. 5 below), and use u a (e) = ijj a {e) in 
(|14|) with 9 a (e) = 0. Then u ap (e) in (|18|) is real and (|19|) coincides 
with (|11|) for resistors. There is no l a = other than l\ = 0, and there 
is no resonance. 
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3 Complex symmetric matrix 

For completeness in this section we give a proof of the theorem which 
asserts (|14|) and determines u a for a complex symmetric matrix. Our 
proof parallels that in ^5] . 

Theorem: 

Let L be an n x n symmetric matrix with generally complex ele- 
ments. Write the eigenvalue equation of \J L as 

L V> a = o~ a ip a , o~ a >0, a = 1, 2, n. (20) 

Then, there exist n orthonormal vectors u a satisfying the relation 

Lu a = X a u* a , a = l,2,...,n (21) 

where * denotes complex conjugation and X a = *Jb T a ~e lda , 6 a = real. 

For nondegenerate a a we can take u a = ip a ; for degenerate a a , the 
u 's are linear combinations of the degenerate ip a . In either case the 
phase factor 9 a of X a is determined by applying \21\) . 

Remarks: 

1. The A Q 's are the eigenvalues of L if real. 

2. If {u a , X a } is a solution of (|21|) . then {u a e lT , X a e 2tT }, r = real, 
is also a solution of ()21j) . 

3. While the procedure of constructing u a in the degenerate case 
appears to be involved, as demonstrated in examples given in section 5 
the orthonormal n's can often be determined quite directly in practice. 

4. If L is real, then as aforementioned it has real eigenvalues and 
eigenvectors and we can take these real eigenvectors to be the u a in 
()21|) with l a real non-negative. 

Proof. 

Since L is Hermitian its nondegenerate eigenvectors ip a can be 
chosen to be orthonormal. For the eigenvector ip a with nondegenerate 
eigenvalue a a , construct a vector 

4> a = (Lip a )* + c a i\) a (22) 

where c a is any complex number. It is readily verified that we have 

L' L 4> a = a a (p a , (23) 
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so 4> a is also an eigenvector of L with the same eigenvalue a a . It fol- 
lows that if a a is nondegenerate then <j) a and ip a must be proportional, 
namely, 

LV« = A Q C (24) 

for some X a . The substitution of (|24j) into p3|) with 4> a given by (|22|) 
now yields |A Q | 2 = a a or X a = ^/a^e iea . Thus, for nondegenerate a a 
we simply choose u a = tp a and use (|2T|) and (|2fl|) to determine the 
phase factor a . This establishes the theorem for non-degenerate X a . 

For degenerate eigenvalues of L, say, a\ = 02 = cr with linearly 
independent eigenvectors ipi and ip2, we construct 

Vl = (LVi)* + V^e iei V'i 

^2 = (LV2)* + \/^e iea ^2 (25) 

where the choice of the real phase factors 0\ , 62 is at our disposal. We 
choose 0i, 62 to make v\ and V2 linearly independent to satisfy 

e^~ 8 ^ = {v2iv 1 )*/{v2,v 1 ) (26) 

where (y, z) = (y T )*z is the inner product of vectors y and z. 
Now one has 

L Vl = V^e^vl, Lv 2 = V^e ie2 v*2 
I/Lui = avi, L V2 = <J V2 ■ (27) 

Write 

ui = «l/H, (28) 

where \v\ = (v, v) is the norm of v, and construct y = V2 — (V2, U\)ui 
which is orthogonal to u\. Next write 

U2 = y/\y\- (29) 

Then, it can be verified by using (|26|) that u\ and 112 are orthonormal 
and satisfy 

Lu 2 = Vae i92 u* 2 . (30) 

In addition, both u\ and U2 are eigenvectors of L with the same 
eigenvalue a, hence are orthogonal to ip a , a > 3. This establishes the 
theorem. 



<s 



In the case of multi-degeneracy, a similar analysis can be carried 
out by starting from a set of v a to construct u Q 's by using, say, the 
Gram-Schmidt orthonormalization procedure. For details we refer to 



If there exist eigenvalues X a = 0, a > 2, a situation which can occur 
at specific frequencies oj in an AC circuit, then the effective impedance 
()19|) between any two nodes diverge and the network is in resonance. 

In an AC circuit resonances occur when the impedances are pure 
reactances (capacitances or inductances). The simplest example of a 
resonance is a circuit containing two nodes connecting an inductance 
L and capacitance C in parallel. It is well-known that this LC circuit 
is resonant with an external AC source at the frequency oj = l/V LC. 
This is most simply seen by noting that the two nodes are connected 
by an admittance y\2 = jojC + 1/jojL = j(ojC — 1/ojL), and hence 
Z12 = I/1/12 diverges at oj = l/V LC. 

Alternately, using our formulation, the Laplacian matrix is 



so that L*L has eigenvalues o\ = 0, 02 = 4| 7/12 1 2 and we have Ai = 
as expected. In addition, we also have A2 = when yi2 = at the 
frequency oj = 1/y 'LC '. This is the occurrence of a resonance. 

An extension of this consideration to reactances in a ring is 
discussed in Example 2 in the next section. 



Examples of applications of the formulation (|19|) is given in this sec- 
tion. 

Example 1. A numerical example. 

It is instructive to work out a numerical example as an illustration. 



D3I- 



4 Resonances 




(31) 



5 Examples 
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Consider three impedances Z\2 = iV%, ^23 = — iV^, 231 = 1 con- 
nected in a ring as shown in Fig. 1 where i = j = y—l. We have the 
Laplacian _ _ 

'l-*/V3 i/VS -1 \ 

i/y/3 -i/V3 ■ (32) 

-1 l + i/VZ/ 

Substituting L into (|12|) we find the following nondegenerate eigenval- 
ues and orthonormal eigenvectors of L, 



01 



(72 



<73 



0, 1/n 
3 - 2\/2, </>2 

3 + 2V2, -03 



24-6\/2 
1 

24 + 6^ 



2 - \/2 + i\/3 
-V2-1 - is/3 

2V2-1 
'2 + V2 + iV3" 
V2-1- iV3 
-2V2- 1 



(33) 



Since the eigenvalues are nondegenerate according to the theorem we 
take Ui = ipi, % = 1, 2, 3. Using these expressions we obtain from (|21j) 



= v^-l, 



■M2 



1 



V2 + I, e ie3 



-i6> 2 



3\/2 - 2 + z\/3(2\/2 + 1) 



3\/2 + 2 + n/3(2\/2 - 1)1 . (34) 



'02 

= 

Now (Unj reads 

^p? = — f^2 P - "29) + ^-7= ( 
^/^2 v v v^3 v 

using which one obtains the impedances 

Z12 = 3 + iV3, Z23 = 3 - iV3~, Z31 = 0. 



-w 3 



. U3p ~U3q , 



(35) 



(36) 



These values agree with results of direct calculation using the Ohm's 
law. 

Example 2. Resonance in a one- dimensional ring of N reactances. 

Consider N reactances jx\, jx2, jxN connected in a ring as 
shown in Fig. 2, where x = ojL for inductance L and x = —1/uC 
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for capacitance C at AC frequency uj. The Laplacian assumes the 
form 



/ vi + vn -yx 
-yx yi + y2 



V 



o 

-UN 










-2/2 
















-yN \ 





-VN-i vn-x + yN -yN 
-yx VN + yx ) 



,(37) 



where yi = The Laplacian L has one zero eigenvalue Ai = as 

aforementioned. The product of the other N — 1 eigenvalues A Q of L is 
known from graph theory |14[I15| to be equal to N times its spanning 
tree generating function with edge weights 2/1,2/2, • ••,2M r - Now the N 
spanning trees are easily written down and as a result we obtain 



N 



1 1 
— + — + 

2/1 2/2 



1 



Ni-j)*' 1 ^ + x 2 + • • • + x N )lx x x 2 ■■■x N . (38) 



H 2/12/2 ■■•VN 

yN, 



It follows that there exists another zero eigenvalue, and hence a res- 
onance, if x\ + x 2 + • • • + 2; at = 0. This determines the resonance 
frequency oj. 

Example 3. A one- dimensional ring of N equal impedances. 

In this example we consider TV equal impedances z connected in a 
ring. We have 

L 

where y = 1/z and 

/ 2 



11 T per 

y n , 



L t = y* T P« 



\y\ 2 (T p N er ) 2 



(39) 



rppcr 



1 



1 



V 


















■1 \ 



-1 

2 J 



(40) 



Thus L and L all have the same eigenvectors. The eigenvalues and 
orthonormal eigenvectors of T^ r are 



Mr! 



2[1 - cos(2nvr/iV)] = 4cos 2 (ravr/iV) 
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1 



2n 



71 



0,l,...,iV-l 



(41) 



(42) 



where u = e l2n / N . The eigenvalues of L are 

<>n = \y\ 2 A- 

Since 

(JN-n = <?n , (43) 

the corresponding eigenvectors are degenerate and we need to con- 
struct vectors u n \ and u r a for < n < N/2. For N = even, however, 
the eigenvalue cr N / 2 i s non-degenerate and needs to be considered sep- 
arately. 

For < n < N/2 the degenerate eigenvectors 

i> n and ^Ar_ n = y>* (44) 
are not orthonormal. Then we construct linear combinations 



U n 2 



V2 



\ 



COS 
COS 



2nir 
N 

N 



, 2(N-l)nn 



y/2i 



sm 
sin 



2nir 

Arm 
N 



, ■ 2(N-l)nn 



n= 1,2,. 



iV - 1 



(45) 

which are orthonormal, where [x] = the integral part of x. The it's are 
eigenvectors of L with the same eigenvalue a n = |y| 2 /Un- F° r N = 
even we have an additional non-degenerate eigenvector 

-1 



u 



N/2 



1 
-1 

V-i7 



(46) 
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We next use (|2*T|) to determine the phase factors Q n \ and 9 n 2- Com- 
paring the eigenvalue equation 

L-Unl = (yHn)u n l with L U n i = (\y\ fJ-n) e^" 1 U* l5 

L"Un2 = {y^n)u n 2 with L u n2 = (\y\p, n )e ldn2 n* 2 , 
and Luat/2 = Kv) u N/2 with = 4|?/|e i6,JV / 2 n^ /2 (47) 

we obtain 



7nl — Pn2 



'JV/2 



(48) 



where is given by y = \y\ e" 

We now use (|19|) to compute the impedance between nodes p and 
q to obtain 



_A_ J_ 



2npir 2nqir \ ' 

cos — — cos ■ ' 



N 



N 



i 2 \ sin 



2npit 2nqir\' 
sin ■ 



N 



N 



E 



(49) 



where [x] denotes the integral part of x and 



E 



1 



2Ny 
0, 



-IF 



, iV = even 

JV = odd. 

After some manipulation it is reduced to 

N—l \gi2npir/N gi2nqn/N 



E 



m N ^ 2[1 - cos(2n7r/iV)] 
This expression has been evaluated in jS] with the result 



Z pq = z\p-q\ 



1 



IP 



N 



(50) 



(51) 



(52) 



which is the expected impedance of two impedances \p — q\z and 
(N — \p — q\)z connected in parallel as in a ring. This completes 
the evaluation of Z 



pq- 



Example 4- Networks of inductances and capacitances. 

As an example of networks of inductances and capacitances, we 
consider anMxJV array of nodes forming a rectangular net with free 
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boundaries as shown in Fig. 3. The nodes are connected by capaci- 
tances C in the M directions and inductances L in the N direction. 



The Laplacian of the network is 

L = (juC) T% cc ® I N - (-L) I M 



where T^ e is the M x M matrix 



i free 
■M 



( 1 





V o 



-1 
















-1 










T free 



\ 





(53) 



2 -1 
-1 1 J 



(54) 



and I at is the N x N identity matrix. This gives 
L*L = ( U ;C) 2 U| f cc ®I J v-2(y) T|, ee ®T^ e +(^-) 2 lM®U^ c (55) 



where U^ e is the M x M matrix 

/ 2 -3 1 

-3 6-4 1 
1-4 6-4 

1-4 6-4 

T Tircc _ 

U M - 1-4 6 




\ 



















o \ 







(56) 



Now has eigenvalues 



-4 6 -3 
1-32/ 



mir 



A m = 2(l-cos0 m ) = 4sin J (0 m /2), 6 m = — 



M 



(57) 



and eigenvector ifj^f ^ whose components are 

m = 0, for all x, 



AM) 
Tmx 



1^ 

M ' 



COS 



(x+^)9 m , m = 1,2,... ,M- 1, for all x. 



(58) 
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It follows that L*L has eigenvectors 
and eigenvalues 

-\r( n ■ 2®rn 1 . 2 0n 

a mn = W\uC sin — - — sin — 
where 9 m = rrnr/M, <p n = nir/N. This gives 



4jLCsin 2 (0 m /2) - _ s in 2 ((/» n /2) 



(59) 



(60) 



an = ±7r/2 • 



(61) 



Since the vectors W^nYCx y) are orthonormal and non-degenerate, ac- 
cording to the Theorem we can use these vectors in (j!9|) to obtain the 
impedance between nodes {x\,y\) and (x2,2/2)- This gives 



Z, 



free 



M-l N-l 
m =0 n=0( mjn )^(o,0) 



\^(m,n);(xi,j/i) r \rn,n);(x 2 ,y2) 



) 



-j . j^L 2j 

\xi - x 2 \ + -^-r \y\ - 2/2 + 



iVwC" 1 M 1 1 MiV 

^l^ 1 ^ 1 [cos (xi + l^9 m cos [yi + - cos (x 2 + |)0m cos (y 2 



m=l n=l 



-uC(l - cos6» m ) + ^(1 - COS0 n ) 



(62) 



As discussed in section 4, resonances occur at AC frequencies de- 
termined from l mn = 0. Thus, there are (M — 1)(-/V — 1) distinct 
resonance frequencies given by 



OJr. 



sin(nvr/27V) 



sin(mvr/2M) 



1 



'LC 



m = 1,...,M - l;n= 1,...,N - 1. 



(63) 

A similar result can be found for anMxJV net with toroidal boundary 
conditions. However, due to the degeneracy of eigenvalues, in that case 
there are [(M + l)/2] [(A r + l)/2] distinct resonance frequencies, where 
[x] is the integral part of x. It is of pertinent interest to note that a 
network can become resonant at a spectrum of distinct frequencies, 
and these resonances occur in the effective impedances between any 
two nodes. 
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In the limit of M, N — ► oo, Q63JI becomes continuous indicating 
that the network is resonant at all frequencies. This is verified by 
replacing the summations by integrals in Q62|) to yield the effective 
impedance between two nodes (xi,yi) and (0:2,2/2) 



(si,3/i);(a:2,3/2) 



4vr 2 



2- 



2tt 



1 - cos[(xi - x 2 )0] cos [(2/1 - 2/2) 
-wC(l - cos6>) + ^(1 - cose?! 



(64) 



which diverges logarithmically/ 



6 Summary 

We have presented a formulation of impedance networks which permits 
the evaluation of the effective impedance between arbitrary two nodes. 
The resulting expression is (|19|) where u a and l a are those given in 
((TU). In the case of reactance networks, our analysis indicates that 
resonances occur at AC frequencies u determined by the vanishing of 
l a . This curious result suggests the possibility of practical applications 
of our formulation to resonant circuits. 
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Fig. 1 An example of three impedances in a ring. 
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Fig. 2 A ring of N reactances. 
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Fig. 3 A 6 x 4 network of capacitances C and inductances L. 
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